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^ Abstract 

^ Two-dimensional conformal field theories with a large central charge and 

r- 1 

a small number of low-dimension operators are studied using the conformal 
block expansion. A universal formula is derived for the Renyi entropies of iV 

> 

disjoint intervals in the ground state, valid to all orders in a series expansion. 
This is possible because the full perturbative answer in this regime comes from 
the exchange of the stress tensor and other descendants of the vacuum state. 

m 

Therefore, the Renyi entropy is related to the Virasoro vacuum block at large 
— , central charge. The entanglement entropy, computed from the Renyi entropy 

by an analytic continuation, decouples into a sum of single-interval entangle- 
ments. This field theory result agrees with the Ryu-Takayanagi formula for 
the holographic entanglement entropy of a 2d CFT, applied to any number of 
intervals, and thus can be interpreted as a microscopic calculation of the area 
of minimal surfaces in 3d gravity. 
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1. Introduction 

Entanglement entropy measures the entropy of a subsystem after tracing out the 
environment. In two spacetime dimensions at a critical point, when the subsystem is 
a single interval, the ground state entanglement entropy is universal in the sense that 
it does not depend on the details of the conformal field theory [U 121 E]- Conformal 
invariance fixes S in terms of the central charge c and the length L of the subsystem. 

The entanglement entropy of a subsystem that consists of multiple disjoint inter- 
vals depends, in general, on the full operator content of the theory [H E]- One aim of 
this paper is to show that in the limit of large central charge, in a theory with a small 
number of light operators, the entanglement entropy is again universal. To leading 
order in 1/c, it is uniquely fixed by conformal invariance, though in a more elaborate 
way than the entanglement of a single interval. Renyi entropies, defined below, are 
similarly universal. 

The motivation for studying this class of CFTs comes from the AdS/CFT cor- 
respondence, since these are the CFTs which are expected to have a holographic 
dual with a good semiclassical gravity description. Entanglement entropy is typically 
difficult to compute in field theory, but in a CFT with a holographic dual, there is 
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a simple and universal proposal [5J. For a single interval in two dimensions, it is 
straightforward to check explicitly that the holographic formula agrees with the stan- 
dard CFT result. In higher dimensions or with multiple intervals, the simple formula 
is a surprise, but it has passed a number of nontrivial tests [7j. A partial derivation 
exists for two intervals in two dimensions, for the first few orders in a series expansion 
[8]. We extend this derivation to all orders and to any number of intervals. 

The result also applies to problems related by a conformal mapping, including 
the (possibly time dependent) entanglement entropy of a CFT at finite temperature 
where it is very natural to consider multiple intervals (9j QUI ITT] . 

The argument relies on a formula for the Virasoro conformal block at large c. This 
formula is well known, but appears to have found few (if any) direct applications in 
AdS/CFT. It is likely that other universal features of 3d gravity, or of 2d CFTs with 
gravity duals, can be understood in a similar way. 

The strategy to compute the entanglement entropy and the relevance of the con- 
formal block are as follows. Divide a system, always taken to be in its groundstate, 
into two parts A and B. The reduced density matrix of region A is obtained by 
tracing out B, pa = TiBptot- This creates a mixed state, with entanglement entropy 

Sa = -Trp^ log p A ■ (1.1) 

One approach to compute this is the replica method. We define the Renyi entropies 

fii B) = -j^ log TipS, (1.2) 

for integer n > 2, then analytically continue n — > 1 to find the entanglement entropy 
Sa = Sa ■ When A consists of iV disjoint intervals, the Renyi entropy can be realized 
as a 2iV-point correlation function. We will expand this correlation function in con- 
formal blocks, and show that the leading contribution to the Renyi entropy at large c 
is captured entirely by the Virasoro block for the vacuum state. Other contributions 
are exponentially suppressed. This statement is true to all orders in the OPE series 
expansion, but fails non-perturbatively as different terms in the expansion exchange 
dominance at large c. Note that unlike higher dimensions, the vacuum block in 2d 
CFT is nontrivial, since it includes the stress tensor and an infinite number of other 
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Virasoro descendants. 

The vacuum block is not known in closed form but can be computed numerically 
by a simple recursion relation to find the Renyi entropy to high accuracy. It can also 
be computed by solving a certain monodromy problem for a second order differential 
equation. In the limit n — > 1, the monodromy problem is solved analytically to 
compute the entanglement entropy. The result for N disjoint intervals is 

S*=lT,**(^) (1-3) 

where z^ for % — 1 . . . 2N is the endpoint of an interval, e is a UV cutoff, and the sum 
is over pairs dictated by the OPE channel in a way described below. This is 
equivalent to the holographic formula [6]. 



The CFT derivation of (1.3) is valid only within some finite region around the 
origin of the OPE in any channel. A complete derivation of the holographic formula 
would require a non-perturbative argument that there are no other phases in the 
parameter space of Zi, i.e., phases which do not include the origin of any OPE channel. 
Of course, the full OPE at finite central charge is convergent and any can be 
expanded in some channel; however, the leading-c term is not analytic in the same 
range of zi as the original OPE. We show that crossing symmetry rules out any other 
phases for N = n = 2 but do not address the general case non-perturbatively. 

The Renyi entropies also match precisely with a gravity calculation done in 



In fact, exactly the same monodromy prescription for computing £7^ that we will 
derive from CFT was derived independently from 3d gravity in a completely different 
way [12J. This means that the semiclassical conformal block with external twist 
operators and an internal unit operator is equal to the on-shell Einstein action on an 
appropriate 3-manifold with nontrivial topology. 

As mentioned above, the technique used to derive the Renyi entropy may have 
other applications in 3d gravity, for example to the calculation of black hole scattering 
amplitudes from CFT. It is a manifestation of the known connection between AdS3 
gravity, SL(2, R) Chern-Simons theory, and Liouville CFT in the classical limit [13l 
HU [15j [16] . We will not explore this triangle of connections in any detail, but return 
to this point of view in the discussion section. 
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2. Semiclassical conformal blocks 



We begin with a discussion of correlation functions in a 2d CFT with large central 
charge using the OPE. The main conclusion will be that in this limit, the stress tensor 
and other Virasoro descendants of the vacuum give a nontrivial contribution to the 
4-point function that can be computed by imposing a trivial monodromy condition 
on a certain differential equation. 

2.1 General operators 

A four-point function of primary operators on the plane may be expanded in conformal 
blocks, 

(O 1 (0)O 2 (x)O 3 (l)O4(oo)) = K h u x)T{c, h p , hi, x) . (2.1) 

p 

We have expanded in the s-channel x — > 0. Here (hi, hi) are the dimensions of Oj=i...4, 
the sum is over primary operators O p of dimension (h p , h p ), and 

a P = Cj2^34: (2-2) 

where c?- is an OPE coefficient. The Virasoro blocks J 7 capture the contribution of 
all the Virasoro descendants. They are unknown in general, but straightforward to 
compute in a series expansion. This can be done very efficiently using a recursion 
formula [T7j described in appendix [Aj 



We are interested in correlation functions at large central charge. In the 'semi- 
classical' limit, defined by taking c large with h{/c and h p /c held fixed, the block 
exponentiates [TSl EH] , 



J-"(c, h p , hi,x) ~ exp 



c , / h p hi 
b \ c c 



(2.3) 



The function / is once again unknown, except in an expansion around x = 0: 
o 2h n 
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However / is determined by the solution of a certain monodromy problem. Consider 
the differential equation 

^"(z)+T(z)i>{z) = , (2.5) 



with 



where (z±, z 2 , z 3 , Z4) = (0,x, 1, 00). The Cj are called accessory parameters. Three of 
them are fixed by the requiring T(z) to vanish as z~ 4 at infinity. This imposes 

y; a = , YVcjZi - — ) = o , ^( C ^ 2 - ^-^o = o , (2.7) 



so that 

c h h 2 h 3 h 1 + h 2 + h 3 - h 4 c c 2 x(l - x) 

6 ~ z 2 (z-xf (l-zf z(l-z) 6 z(z - x)(l - z) ' 1 J 



The differential equation (2.5) has two solutions, ipi and ip 2 . As we take the solu- 
tions on a closed contour around one or more singular points, they undergo some 
monodromy 

O—O- (2 - 9) 

The 2x2 monodromy matrix M depends on the basis of solutions ipi t2 , but its trace 
is invariant. 

The connection between this differential equation and the semiclassical conformal 
block is the following. First, we choose the accessory parameter c 2 (x) so that the 
monodromy on a cycle enclosing both and x is 

TrM 0:c = -2costtA p , h p = ^(1 - Aj) . (2.10) 

Then the semiclassical Virasoro block is determined by 

g = *M. P.ll) 
The integration constant in this equation is fixed by comparing to the series expansion 



Although the formula (2.3) has never been proved directly from the definition 
of the Virasoro block as a sum over descendants, it follows from the path integral 
of Liouville theory, and has passed extensive checks. We will review the Liouville 



derivation in section 4.1 below and generalize it to more than four external operators. 

The Virasoro conformal block depends only on the algebra, so this result from 
Liouville theory is applicable to any CFT in the semiclassical limit. 

2.2 Unit operator 

The semiclassical limit c — > oo is usually taken with h p /c held fixed. In fact, if the 
external weights are equal, hi = h, then we can also take the semiclassical limit with 
hp held fixed, i.e. 7 = h p /c — > 0. The limits commute: 

lim lim - log J-"(c, 07, c(3, x) = lim - log J-"(c, 0, c(3,x) . (2-12) 

7— >0 c— »oo c c— >oo c 

This can be seen from the recursion representation of J 7 in appendix [Af*] Let us denote 
the semiclassical block for these light operators by 

fo(-,x] =f(o,-,x] . (2.13) 



C J \ c 

The dependence on h p disappears, so we will set h p = and refer to this as the 
vacuum block. In d > 2 spacetime dimensions, the vacuum block is trivial; it gives 
only the disconnected piece of a correlator. However in d = 2, the vacuum block 
includes the stress tensor and all of the other Virasoro descendants. 

It can be computed (in principle, or numerically) by solving the monodromy prob- 
lem above, where now we tune C2(x) to impose trivial monodromy around the points 
0,x: 

TrM 0x = 2 . (2.14) 



This follows from (2.10) with h p = 0. The semiclassical vacuum block can also be 
interpreted as a classical on-shell Liouville action; we will return to this below. 

*In the recursion, the series expansion of log J 7 is organized so that the only h p dependence is a 



sum over terms of the form ^ ° b ^ fc with b — O(c). The limits in (2.12) commute term by term. It 
is enough for the external weights to be equal in pairs, hi = h 2 , h 3 = /14. 
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3. Entanglement entropy of two intervals 

3. 1 CFT calculation 

Suppose that we have a family of theories, labeled by the central charge, that admits a 
large-c limit. Besides unitarity and compactness, we will make only two assumptions 
about this family of CFTs. First, in the large-c limit, correlation functions are smooth 
in a neighborhood of coincident points and obey cluster decomposition. This restricts 
the growth of OPE coefficients to be at most exponential in c; otherwise the vacuum 
does not appear in the OPE 0. Second, we assume that the number of light 

operators does not grow with c. Specifically, the density of states is d(A, A) = O(c ) 
for A, A < A gap where we will take A 5ap = c/24 but this can be relaxed to A gap = 
0(c). Examples are symmetric product theories and theories with a gravity dual (in 
this case A gap is the energy of the lightest black hole). We will compute the Renyi 
entropies and entanglement entropy in theories of this type to leading order in 1/c, 
without using the specific details of the CFT. 

The replica method, described in the introduction, is a useful approach to entan- 
glement entropy when the system is in a state where the wavefunction is computed 
by a path integral on some manifold M. This includes the vacuum state, or a thermal 
state in CFT. In this case the replica partition function can be characterized in two 
equivalent ways. First, it is a path integral on an n-sheeted cover of M, with branch 
cuts on region A where the sheets are connected. This manifold is a singular Riemann 
surface with nontrivial topology. Second, the partition function can be computed in 
an n-fold product theory on the original manifold M, with twist fields inserted at the 
boundaries of region A. The twist fields are defined to glue together the n copies of 
the CFT in a way that reproduces the path integral on the higher genus surface (see 
[3] for a review) JT] 

^It is not guaranteed that the replica trick will always produce the correct entanglement entropy 
because of ambiguities in the analytic continuation n — > 1. We will assume without proof that it 
works with the obvious choice of analytic continuation in the CFTs we consider, and furthermore 
that we can perform the analytic continuation on the leading-c term directly, but we caution that 
this is not entirely justified. 
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For region A consisting of N disconnected intervals, 



A=[z 1 , z 2 ] U [z 3 , z*] U • • • U [zzn-i, z n] 



Z\ < Zi < ■ ■ ■ < Z2N 



(3.1; 



the Renyi entropy, computed in the product theory on M, is [21 [3] 



exp((l-n)5? ) ) = ($ + (^ 1 )$_(^ 2 )...$ + (^ iV _ 1 )$_(^ Ar )) . 



(3.2) 



The twist operators Q± have weight 



(L ,L ) = (H,H) , H 



24 



n 



1/n) . 



(3.3) 



For a single interval, S 1 ^ is completely fixed by conformal invariance, but for N > 2 
it generically depends on the full operator content of the theory. 

In the rest of this section we set N = 2. This case has been discussed in detail in 
[8]. We will repeat parts of that discussion in our language, both in order to setup the 
problem of general N and to clarify and extend some aspects of the CFT calculation 
in 0. 

We can set (z\, z%, z%, z\) = (0, x, 1, oo) by a conformal transformation. The cross- 
ratio x is real. The Renyi entropy is a 4-point function, so it has the conformal block 



expansion (2.1). At large central charge, expanding in the s-channel, 



exp ^(1 — n)S^^j = a p exp 



o \nc nc 



nc 
~6 



f 



h p H 

nc nc 



x 



, (3.4) 



where we have used the central charge of the replica theory, nc. We want to show that 
the exponential dependence is entirely captured by the first term in the sum; low-lying 
terms affect only the irrelevant coefficient of the exponential, and high-dimension 
terms are non-perturbative in 1/c. This is typical behavior of a thermodynamic 
function with different phases. The argument is simplest if we assume a p does not grow 
exponentially with c^\ Then our assumption about the low-lying operator spectrum 



allows us to ignore the coefficient and any multiplicities in (3.4) for h p , h p < A gap . In 



■^This extra assumption is not necessary, though without it the region in which the phase (3.7) 



dominates may cover a smaller range of x. The assumption that we have smooth correlators in a 
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this range, f(S p , —,x) is an increasing function of 5 P for any x < |, so the low-lying 



terms exponentially dominate. Similarly, heavy operators in (3.4) are suppressed 
non-perturbatively in 1/c. Therefore in a neighborhood of x = 



S 



(n) 



nc 



3(n 



tv/o ( — ,x 
1) \nc 



(3.7) 



This is one of our main results. It is the contribution from all light operators, or 
equivalently the contribution from the vacuum and its descendants. We reiterate 



that (3.7) is the full answer for the Renyi entropy to leading order in 1/c in a finite 
region around x = 0, i.e., to all orders in a series expansion in x. We will return 
to the question of how big this region is below. There are both perturbative and 
non-perturbative corrections in 1/c. We have not used any special properties of the 
twist operators other than the fact that H ~ c/24, so a similar formula applies to the 
leading-c correlation function of other heavy operators. 



The result (3.7) allows for a simple numerical or series calculation of the Renyi 
entropy for any n, not necessarily an integer. This can be done easily using the 
recursion formula in appendix [A] to compute log J 7 . The results for various n are 
plotted in figure [l|a), and the series expansion is given in the appendix. Alternatively, 
it can be computed by solving the trivial-monodromy problem that defines /o, but 
for general n this is more difficult than just computing log J 7 directly. 

The advantage of the monodromy prescription is that it can be solved analytically 
in the limit n — > 1,H — > relevant for the entanglement entropy. Denoting n = 



neighborhood of x — as c — > oo already implies a bound of the form 

a p < exp 





' hp 


hp\~ 


I 9 


\nc 


nc J 



(3.5) 



Therefore, up to an irrelevant pref actor to account for multiplicities, we can bound the sum over 



low-dimension operators in (3.4 1 by 

-1/24 /-1/24 



dSr, 



dS p exp 



\ (g($p, 8p) - f(S p , —,x)- f(S p , —,x) 
D \ nc nc 



(3.6) 



This is exponentially dominated by an endpoint of the limits of integration, or by a saddlepoint; 



near x = 0, it must be dominated by 8 P 



0, because the vacuum must give the leading term 



in the correlation function near coincident points. 
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Or 





Figure 1: The vacuum contribution to the Renyi entropy (divided by c/6), for 
n = 2,5, 10, 40, from bottom to top. In (a) we show the contribution to the s-channel, 
and in (b) both channels. These are computed from (3.7) using the recursion relation 



in appendix [X] iterated 10 times. The error, estimated by adjusting c and discarding 
the last term in the recursion, is less than 1% for z < 0.99. The loge UV divergence is 
dropped. 



1 + 12a, the entanglement entropy is 



c 1 
Sa= ^ lim -/oCa, a) 

36 «->0 a 



(3.8) 



so we only need to solve the monodromy problem with vanishingly small operator 



insertions. For small hi in (2.6), the accessory parameter C2 will also be small, so 



all terms in T(z) are important only at the singular points. Therefore to solve the 
monodromy problem around 0, x, we can ignore the other singular points, and consider 
the simpler problem with 



T(z) = 6a — + 



(z — x) : 



+ 2 ) + 



X 



(-2 



z — X 



(3.9) 



This has trivial monodromy around 0, x if T is regular at infinity. Therefore we 
require the sum of the residues of simple poles to vanish: 



C-2 



12a 



x 



+ 0(a 2 



(3.10) 



Integrating, 



fo(a, x) = 12a log x + 0(a z 



(3.11) 



This is the leading-c contribution to the Virasoro block for a heavy internal operator, 
analytically continued to small h v \ note that it is not equal to the block for a light 
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internal operator. The formula (3.11 ) states that only the identity operator - not the 



stress tensor or its descendants - contributes, at any value of x. Thus the entangle- 
ment entropy for two intervals will factorize, although this was not true of the Renyi 



entropies with n > 1. From (3.8) we find the entanglement entropy 



c fx 
3 l0g U 



(s-channel) 



(3.12) 



(In this expression and similar expressions below we have reintroduced the UV cutoff 
e that is necessary to regulate the twist operators.) 

We chose to expand in the s-channel x ~ 0, but now let us expand in the t-channel 
x ~ 1. Blocks in the t-channel are related to those in the s-channel by x — > 1 — x 



so the answer can be obtained from (3.12) but let us demonstrate this directly. The 



argument is identical, except that now we define fo by imposing trivial monodromy 



around a cycle enclosing x, 1 instead of 0, x. Instead of (3.9) we have 



T(z) = 6a ( - — + ' 



2 \ r-j.r r 2 



[z — l) 2 [z — x) 2 z — 1 / z — 1 z — X 
Regularity implies c 2 = 12a /(x — 1), and dfo/dx = c 2 gives 



(3.13) 



fo = 12cdog(l -x) + 0(a 2 



(3.14) 



Therefore 



(i-channel) , 



to all orders in a series expansion around x = 1. 



(3.15) 



The results (3.12, 3.15) were derived to order x by a series expansion of the 



Virasoro blocks in [8j. From the present derivation, which follows similar logic but 
exploits the formula for the semiclassical blocks, it is valid in finite regions around 
x = and x — 1. 



However we have not specified the range of x in which (3.7, 3.12, 3.15) are valid. 



To do this we would need to determine when the heavy operators in (3.4) first start to 



dominate. The s-channel result (3.12) must break down at or before x = |, because at 



this point, the heavy operators in the s-channel must account for the vacuum running 
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in the t-channel. The question is whether there is another transition at some point 
x = x c , with < x c < |, to a phase that does not correspond to the vacuum block in 
either channel. (Other phases certainly exist elsewhere on the complex plane, but in 
our case x is real and < x < 1.) Such a phase would not be visible in a perturbative 
expansion (in 1/c) around the origin of the OPE. 

The answer for n = 2 is that the s-channel Renyi entropy is valid for < x < |, 
and there is a sharp transition at x = | to the t-channel result. There are no other 
phases. This was demonstrated in [5] by mapping the n = 2 twist correlator to 
a torus partition function. An alternate derivation is given in appendix [B] using 
crossing symmetry of the OPE to show that the contribution of heavy operators in 
the s-channel is equal to the vacuum term in the t-channel. The argument is limited 
to n = 2 unless we add the assumption that certain OPE coefficients are exponentially 
suppressed at large c, so we leave open the question of whether there are other phases 
in the higher Renyi entropies or the entanglement entropy. It may depend on the 
particular theory. If we assume that the dominant contribution is always the vacuum 
in some channel, then the entanglement entropy is 

Sa = jjmin jlog (?) , log H-^j } , (3.16) 

and the Renyi entropies for various n are as shown in figure [T^b). 

3.2 Comparison to holography 

In [B] it was conjectured that the entanglement entropy of a CFT with a holographic 
dual can be computed by a simple geometric formula. The proposal is 

S A = mm 7 — — — , (3.17) 

where 7 is a surface of dimension d — 1 in a fixed-time slice of AdS^+i that meets the 
boundary of region A at the boundary of AdS. In AdS3, Newton's constant is Gn = 7^., 
and the 'area' is the length of a spacelike geodesic. When A is a disconnected region, 
the bulk surface can also be disconnected, and it is the total area that computes the 
entropy. 
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The result when A is two disjoint intervals agrees with the CFT calculation above 
[8]. We have two choices for how to draw the geodesies that end on the endpoints of 
region A: 



or 




The configuration that gives the entropy is the one in which the geodesies have min- 
imal total length. This is the disconnected one for small cross-ratio x and the con- 
nected one for large x, with a sharp transition at x — |. Note that in any particular 
channel, the answer decouples into a sum of single-interval entanglements. From the 
CFT point of view, this reflects the decoupling of the different monodromies that must 
be imposed to compute the semiclassical conformal block, when the singularities of 
the differential equation have vanishing weight as n — > 1. 

The holographic entanglement entropy of a single interval gives the classic CFT 
result mentioned in the introduction, = | log ( 2 '~ Zj j , From the diagram, the 
holographic result for two intervals is 

S A = min (5*12 + 5*34 , S 14 + S 23 )- (3.18) 

This agrees with the CFT calculation, applied to the entire range < x < 1. Assum- 
ing the holographic formula is correct, this suggests that in these CFTs the leading 
contribution to the Renyi entropy indeed always comes from the vacuum block in 
some channel. A closely related conjecture about the classical Einstein action on 
manifolds with a higher-genus boundary was made in [19] . 



The Renyi entropies given in (3.7) also agree with a holographic calculation, 



performed in [T2] . The match provides a very simple geometrical interpretation of the 
semiclassical vacuum block /o with external twist operators: it is the Einstein action 
of a particular 3-manifold. We will review the gravity calculation briefly. See [12] for 
details. 

The Renyi entropy is the partition function on a singular Riemann surface M 2 of 
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genus n — 1. In AdS/CFT, this partition function is related to the on-shell Einstein 
action of a spacetime M3 whose conformal boundary is M%. In general, there are 
many manifolds M3 with the same conformal boundary, and we must choose the 
dominant saddlepoint. Within a finite region around the origin of any OPE channel, 
we may restrict our attention to a certain class of handlebodies: given a genus-g 
Riemann surface, we choose g independent cycles ^i,..., g to 'fill in' to construct the 
bulk spacetime. The on-shell Einstein action on a 3-manifold M 3 constructed in this 
way is equal to an appropriately defined Liouville action on M 2 [2U1 ED E21 E3E], 

Q-SEinstein(Ms) _ g - | S ZT ( M2 ,Tj ) 

The Liouville action Szt, defined in [23], depends on M 2 and on the choice of con- 
tractible cycles. It was demonstrated in [12], using [2H E31 EE], that Szt is the 
generating function of the accessory parameters for the trivial-monodromy problem 



defined in section 2.2 dSzr/dzi = riCi. This is exactly the same prescription de- 
rived from the conformal blocks, so Szt — n fo- It follows that the Renyi entropy 
computed from CFT agrees with the gravitational partition function, with classical 
Liouville providing a natural link between the two calculations]^] The choice of trivial- 
monodromy cycles in the definition of fo, i.e. the choice of OPE channel, corresponds 
to the choice of cycles on M 2 that are filled in to construct the spacetime M 3 . 

4. More than two intervals 

In this section we generalize the discussion to the entanglement entropy of N > 2 



intervals. Our aim is to compute the 2iV-point twist correlator (3.2) at large c, to 
all orders in a series expansion around coincident points. This is conceptually very 
similar to the case N = 2: we find a monodromy prescription for the large-c conformal 
blocks, then analytically continue n — > 1 where the monodromy problem can be solved 
analytically. 

§The semiclassical block / also appears in other contexts, including the Bcthe ansatz and 4d 
gauge theory. See [26j [27] and references therein. 
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4-1 Semiclassical k-point conformal blocks 
First we must show that any fc-point conformal block in the semiclassical limit can 



be computed by solving a monodromy problem similar to the one in section 2.1 We 
will derive this from the path integral of the Liouville CFT at large c, generalizing 
the argument for the 4-point function [T7j as reviewed in [281 EH]. The blocks are 
determined solely by the Virasoro algebra so the result applies to any CFT at large 
c. 

First we describe the final result, since this is the only fact needed for the rest of 
the paper. It is the obvious generalization of section |2~T A general Appoint function 
on the plane 

(O 1 (0)O 2 (z 2 )O 3 (z 3 ) • • • O fc _ 1 (l)O fc (oo)) (4.1) 

can be expanded in many different channels. The channels are represented by tree 
graphs with propagators and 3-point vertices. An example for the 6-point function is 
the channel „ Q 



(4.2) 



O p (0) O,(0) O r (l) 

1 6 5 

where p, q, r label primaries appearing in the OPEs. The conformal block in this 

channel is defined so that 



(4.3) 



p,q,r 



More generally, we draw a tree diagram with k external points marked 1 ... k in 
cyclic order, and label each internal line by a primary operator. The blocks can be 
computed in a series expansion by summing over Virasoro descendants up to a given 
level J3H1 EI] , as described for example in [321 [33] . 

In the semiclassical limit c — > oo with hi/c held fixed, the /c-point Virasoro block 
in any channel exponentiates, 



T ~ exp 



c I h„ h 



(4.4) 
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where hi denote the external weights and h a the internal weights. The semiclassical 



block / can be computed as follows. Consider the differential equation (2.5) with T(z 



as in (2.6 ), now summing over i — 1 . . . k. Regularity imposes three conditions (2.7), so 



there are k — 3 accessory parameters, which can be used to tune the k — 3 independent 
monodromies of the differential equation. Choose cycles around the singular points 
which correspond to the OPE contractions in the chosen channel. That is, for each 
contraction Oa(za)Ob{zb) — > Oc(za), we choose a contour ^ c enclosing z& and z B . 
In the 6-point example above, the cycles can be chosen as follows: 




(4.5) 



Now, choose the accessory parameters Cj = 2,..,,fe-2 so that the monodromies on these 
cycles are 

TrM a = -2cos7rA a , h* = ^(1 - A 2 a ) , (4.6) 

where h a is the weight of the primary operator appearing in the corresponding OPE 
channel. Finally, the semiclassical block is computed by integrating 



(4.7) 



We now restrict to equal external weights hi = h. What is the analog of the 
'vacuum block' that was important in the previous discussion? This depends on the 
channel. In some channels, we can choose every internal weight to vanish, 



fa 



h 



f 



h 



(4.8) 



This is possible only if the external operators are contracted in pairs, such as in the 
channels 



and 




In a channel where an external operator meets two internal operators, we cannot 
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set both internal operators to the identity, since this would vanish. In this case, we 

contract into the identity operator whenever possible, and set other internal weights to 

h. This is what we mean by the 'vacuum block' fo in the general case. It can always 

be computed by imposing trivial monodromies on the appropriate cycles. This is 

illustrated by the 6-point example above, where the vacuum block has weights 

3 



1 










h 






(4.9) 

This is computed by imposing trivial monodromy M = \ixi on the cycles 7 P and j r 



in (4.5) (which automatically fixes the monodromy on 7 g to correspond to h). 

This procedure can be implemented numerically, though we will not need the 
explicit results. 

In the rest of this subsection we will sketch a derivation of these statements, 
closely following the logic used for the 4-point function [T71 EHl EH] . The idea is to 
consider the Liouville correlation function {\0\ ■ ■ ■ Ok), where \ is a light operator 
corresponding to a null state, (x\x) — 0. The null decoupling equation for this 



correlator is the differential equation (2.5), and the monodromy condition projects 



onto the contribution from a particular conformal family. 

In Liouville theory (see [3^ [29j [28] for reviews that we follow closely), it is conve- 
nient to parameterize the central charge as 

c=l + 6Q 2 , Q = b+\. (4.10) 

b 

Local operators V a = e 2a< ^ have dimension A a = a(Q — a). Correlation functions are 
defined by the path integral 



Gfc(zi) = (V ai (*i) ■ • ■ V ak (z k )) = J D<pV ai { Zl ) ■ ■ ■ V ak {z k )e~ Q s ^ (4.11) 
with the Liouville action 

S L [<t>) = i- J d 2 z (\d(j)\ 2 ) + fie 4 ') + boundary terms . (4.12) 
In the semiclassical limit b — > 0, correlation functions of heavy operators are given by 
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the action of a classical saddlepoint, 



Gk(zi) ps e s 



S c i 



(4.13) 



where the saddlepoint obeys boundary conditions at the insertion points that depend 
on hi/c. 

Let us choose an OPE channel, and denote the contribution from a particular 
choice of internal weights by Gk(zi)\ M; _„. That is, 



Gk(z. 



k \"i I \p,q,.- 



(4.14) 



where c a is the OPE coefficient computed exactly by DOZZ [3"5"j I3"4"]. (The other two 
indices on c a are suppressed, and depend on the channel.) In the semiclassical limit, 
Gk on the l.h.s. and the c a on the r.h.s. can all be represented as e~^ Scl , so the blocks 



T have the exponential form (4.4). 



The monodromy prescription for the semiclassical block comes from studying the 
decoupling equation for a light null state. The light operator ijj with = — \ — \b 2 
has a null state at level two, 



Ix) 



2(2A^ + 1) 

If we insert this operator into the fc-point correlator, 



(x\x) = o 



(4.15) 



9=(iKz)V ttl (z 1 )---V ah (z h )) 



(4.16) 



then decoupling of the null state implies [18 



g2 2(2Ay, + 1) ^ 



i=l 



[Z — Zi) 



+ 



[Z — Zi 



* = 



(4.17) 



Expanding the heavy operators in conformal blocks, denote the contribution of a 
particular set of primaries by ^f p>qi .... The different ^ p , q ,... have different monodromies 



as we move z around the Zj, so (4.17) must be true for each contribution individually. 



Now we will evaluate ^ p , q ,... semiclassically to simplify (4.17). The insertion of a 
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light operator does not change the classical saddlepoint, so the path integral gives 



* M ,„. « ^{z)c p c q . . . exp (-|(/(z<) + /&))) , (4.18) 

where ^ is a classical field that can be interpreted as the expectation value of ip in 
the presence of the heavy operators. Plugging into the decoupling equation ( |4.17[ ), it 



becomes the differential equation ( 2.5|2.6 ) with accessory parameters q = df/dz{ as 
claimed. 

The last step is to relate the monodromies of ip{z) to the weights h a of the operators 
running in the OPE. Suppose the channel involves a contraction O a Ob — > O c . Then 
if ip(z) is inserted somewhere on a cycle enclosing O a and Ob that is taken small 
enough to avoid other operators, when we do the OPE to produce the conformal 



block expansion of (4.16) we come across (O a (z a )Ob(zb)ip(z)O c (z c )). This correlator 



satisfies a null decoupling equation like (4.17). Plugging into this decoupling equation 



the leading OPE behavior ip(z)O c (z c ) ~ (z — z c ) K Od(z c ) and looking around z ~ z c 
fixes 



|(1 ± Ac) , where h c = — (1 - A*) . (4.19) 
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Finally, the monodromy of the two solutions (z — 2; c )2^ 1±Ae ^ on a cycle enclosing z c is 



g -inA c 




M c = -\ , TrM c = — 2cos7rA c . (4.20) 



This is the same as the monodromy around z a and Zf,, so this completes the argument. 

4-2 Entanglement entropy 

Given the monodromy prescription for A;-point conformal blocks, the discussion of the 
entanglement entropy for intervals is essentially identical to that for two intervals 



in section pi Region A is taken as in (3.1). The twist correlator (3.2) is expanded in 



conformal blocks in a particular channel with k = 2N external operators. 

The leading contribution in this channel comes from light operators exchanged 
between the coincident $ + and We can replace these operators with the identity, 
because the semiclassical block depends only on h p /c — » and the overall coeffi- 
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cient will be irrelevant when we take the log. Therefore the leading non-vanishing 
contribution is what we defined as the vacuum block in the previous subsection, 

3(n — 1) \nc J 

Because the blocks depend exponentially on c, in the large-c limit this is the full 
answer to all orders in a series expansion in this channel. 



The conclusion is that for any channel, the Renyi entropy is given by (4.21) in a 
finite range of Zi, with fo defined by imposing trivial monodromies in a way dictated 
by the channel. As with iV = 2, we have not ruled out the possibility that there are 
ranges of Zi where the dominant contribution is not the vacuum block in any channel. 



The numerical calculation of the Renyi entropies from (4.21) is straightforward but 
computationally intensive, since we do not have the advantage of a fast recursion 
formula as we did for the four-point blocks. 

The derivation of the entanglement entropy proceeds as in section [3} We must 
impose trivial monodromies on certain cycles of the differential equation ( 2.5|2.6 ) 



with external weights hi = H defined in (3.3) and take n — > 1. For n — > 1, H — > 



all of the terms in T(z) are important only at the singular points, so the problem 
decouples into N independent monodromy problems. In each case, if we pair $+(z,) 
with then setting n = 1 + 12a we find Cj = 12a/ (2^ — Zj) + 0(a 2 ). Therefore 

the semiclassical vacuum block is 

/o(a, zi) = 12a £ logfc - z 3 ) + 0(a 2 ) , (4.22) 
and the entanglement entropy is (restoring the UV cutoff) 

^ = iE lo s(^) • (423) 

The sum is over pairs dictated by the OPE channel, selected as follows. First, 
we pair operators which are directly contracted in the OPE; in the 6-point example, 
these are (1, 2) and (4, 5). Then, we pair operators connected in the OPE tree diagram 



by a string of internal if's, like the pair (3,6) in (4.9). These are the correct pairings 
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(a) (b) 

Figure 2: Example of the entanglement entropy of 5 disjoint intervals at large 
central charge in a particular channel, (a) in CFT and (b) holographically. In 
(a), trivial monodromy is imposed on cycles in the z-plane corresponding to 
the dashed lines. The solid semicircles in (b) are geodesies in AdS3 with radial 
direction r, and A is the shaded region on the boundary. 

because these are the ways we can divide the OPE tree diagram cutting only identity 
operators, so these correspond to the cycles of trivial monodromy. 



The result (4.23) is in a particular channel. Each channel is an upper bound for 
the full answer. If we assume that the regions where the vacuum block dominates in 
various channels cover Zi everywhere on the real line, then the full answer is obtaining 
by selecting the channel where this is a minimum, 

S A = mm C -J2 l °z( Z - 1 ^ 3 -) ■ ( 424 ) 



The result (4.24) is identical to the holographic formula for entanglement entropy 
[5J. Each term is ^k— = | times the length of a geodesic through a constant-time 
slice of AdS3, connecting the two twist operators. 

A more complicated example is illustrated in figure [2j in bulk and boundary. In 
the CFT, we draw the OPE tree diagram, and place a unit operator '0' or twist 
operator 'if' on each internal line. Then we pair the twist operators by drawing 
cycles which are only allowed to pierce internal lines carrying a unit operator. The 
semiclassical vacuum block fo is computed by imposing trivial monodromy on the 
corresponding cycles in the z-plane. Finally, the analytic continuation to n — 1 equals 
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the entanglement entropy computed by the length of geodesies in AdS 3 connecting the 
paired operators. As a consistency check, we confirmed numerically that accessory 



parameters obtained by differentiating (4.22) indeed lead to trivial monodromies in 
this 10-point example. 



5. Discussion 

The universal behavior of entanglement entropy in 2d CFTs with holographic duals 
is in some ways an extension of Cardy's formula for the thermodynamic entropy 
|36j . The Cardy formula holds universally at high temperature T — > oo, but in 
theories with a holographic dual, it applies above a critical temperature T c — 1. The 
sharp transition at T c is a consequence of the small number of light operators. It is 
holographically dual to the Hawking-Page phase transition between black holes and a 
thermal gas in AdS 3 , and the Cardy formula equals the black hole entropy [3TI [38| [39] . 
This gives a statistical origin of black hole entropy in many examples. 

According to the Ryu-Takayanagi formula, entanglement entropy is a way of gen- 
eralizing this success to geometric surfaces other than black hole horizons. Patching 
together spacetime via entanglement may be a useful way to approach quantum grav- 
ity [10], and if so, these minimal surfaces must play an important role. The universal 
formula for the entanglement entropy of a 2d CFT provides the microscopic origin for 
the area of these surfaces. Like the Cardy formula, it is really only semi-microsopic, 
in the sense that it relies on a quantum field theory with a microscopic definition but 
does not involve the detailed description of states in that theory. 

The gravitational sector of the AdS / CFT correspondence in the semiclassical limit 
can be summarized by the statement that gravity in d + 1 dimensions is the theory of 
stress tensors in <i-dimensional CFT. That is, gravity describes the thermodynamic 
sector of the CFT. In d = 2, the stress tensor lives in the same Virasoro representation 
as the vacuum state, so with no other matter present this becomes the statement that 
3d gravity at weak coupling is the theory of the Virasoro vacuum representation at 
large c. On trivial topology this is the fact that gravitons are Virasoro descendants 
|41j ; on a torus it implies the Cardy formula for black hole entropy; and one con- 
sequence at higher genus is (modulo the caveats in the derivation) the holographic 
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formula for entanglement entropy. 

Semiclassical Liouville theory provides a link between these two descriptions, grav- 
ity and CFT. On the one hand, the on-shell action of 3d gravity is equal to a Liouville 
action [T51 HS1 ED] • On the other hand, since the thermodynamics of these CFTs is uni- 
versally determined by the algebra, it can be computed in Liouville theory at large c. 
Note that this does not mean quantum Liouville theory is equivalent to either gravity 
or the CFT (see |I2J|13] for obstacles to this interpretation). Rather, Liouville is use- 
ful because in the semiclassical limit it naturally computes thermodynamic quantities 
on both sides |44j . 

To be more specific, a natural conjecture extending our results is that in the 
class of theories we have considered, the leading-c partition function on any Riemann 
surface M 2 is 

Q 

log Z C ft = -- mm S Z t(M 2} T) , (5.1) 
o r 

where T indicates the choice of cycles necessary to specify the Zograf-Takhtajan Li- 



ouville action discussed briefly in section 3.2 See [15] for related comments about 



universality in these theories. It would be very interesting to prove this statement, 
and to understand exactly what restrictions must be placed on the class of CFTs to 
do so. 

This discussion is of course special to two-dimensional CFT. In d > 2, the vacuum 
representation is trivial. The analogous calculation would be to include the exchange 
of all operators built from the stress tensor and its derivatives. This is a much harder 
problem, because there is no local d-dimensional field theory that captures this sector 
of the CFT, like Liouville does in two dimensions. On the other hand, in a 3d CFT 
with higher spin symmetry, correlation functions of currents are fixed by symmetry 
[46] . This is similar to the situation in 2d CFT, so it is plausible that a Liouville-like 
classical theory connects these theories to their higher-spin gravity duals (see jH] and 
references therein). 

Supplemental material 

A Mathematica notebook implementing the conformal block recursion, the numerical 
calculation of monodromies and the accessory parameter c%, and the consistency check 
of the lOpt example in figure [2] is provided with the arXiv submission for download. 
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A. Calculations of the conformal block 

In this appendix we discuss the calculation of the 4-point Virasoro block J 7 and the 
semiclassical block /, and compute the Renyi entropy for two intervals in a series 
expansion. A Mathematica notebook implementing the recursion relation and the 
numerical monodromy algorithm is provided with the arXiv submission of this paper. 

The Virasoro block may be computed by a brute-force series expansion, but this 
tends to be the slowest method for numerical calculations. Two efficient recursion 
formulae exist, one using the analytic structure as a function of c [IE] and the other 
as a function of the internal weight h p [T7] (see also [31]). The latter converges very 
quickly for hi ~ c, so the first few terms can be used to quickly compute T for, say, 
\x\ < 0.99. It is an expansion in 

q{x) = e -*m-*)/K(x) (A.l) 
where K is the complete elliptic integral of the first kind. The inverse is 

(Note that this is the same as the transformation to the torus in [I9j[8]). The recursive 
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formula for the 4-point conformal block on the sphere is: 



l-c c-1 



T{c, h p , hi, x) = (16q) hp+ 24 x 



-h 1 -h,2i 



1 — x) 24 



-h 2 -h 3 



eM^'^H^h^Kq) 



H(c,h p ,hi,q) = 1 - \ ^ IT— FT — "t-T" \ H ( c i H mn + ran, h h </) 



m,n>0 



B m n ( hp H n 



a±{c) 



24 



l-c /25-c 



24 



(A.3) 



c-1 
24 



+ 5 



2 



A mn (c, /ij) — J [ J^J(^3 + A4 — 5'p g )(A3 — A4 — Spg) 

X (Ai + A2 — S pq )(\2 — Ai — Spg) 

5 m „( c )= n n 2s ab 



V 1 



m n 



a=— m+l b=— n+1 



Xi(hi,c) = \ hi + 



1-c 
24 



S rs (c) = + as) 



Zeroes are omitted in the product B mn . The product in A mn is taken over p = 
—m + 1, — m + 3, . . . , m — 1 and q = —n + 1, — n + 3, . . . , n — 1. For h p = with 
equal external weights, (m,n) = (1, 1) is omitted in the sum on the second line. For 
analytic calculations, it is useful to parameterize c = 1 + 6(6 + 1/b) 2 . 



The monodromy problem described in section 2.1 can also be solved numerically to 
compute the semiclassical block / or / . This is done by computing the monodromy 
TrMoz at fixed accessory parameter ci, solving for C2 by shooting, then integrating 



with the boundary condition (2.4) to obtain /. For a fixed value of the accessory 



parameter c-i, the monodromy invariant of the differential equation can be computed 
by numerically solving the ODE on a circle in the z-plane with two arbitrary initial 



conditions to find ipi t 2- Equation (2.9) and its derivative fix the entries of M. 



Now we will use the recursion formula to compute the Renyi entropy of two in- 
tervals in a series expansion around x ~ at large c. With central charge Cq, equal 



external weights hi = cq5 and internal weight h p = 0, from (A.3) we find J 7 



e e 



25 



with 



fo(S,x) = 125 log x + a 2 q 2 + a 4 q 4 + a 6 q 6 + 0{q 8 ) (A.4) 
a 2 = -30725 2 

5 

4096 

a 6 = 5 2 (195 - 184325 + 7843845 2 - 145489925 3 + 985661445 4 ) 

35 

The Renyi entropy is 

(„) = en / 1_ _ 1 \ 
A 3(n- 1) V 24 n 2 / 

For comparison to [8], the Renyi mutual information is defined by subtracting the 



log term from (A. 5). The result is equation (4.41) of [S]; when n = 2 it agrees with 



equation (4.30) of [S]. 

B. Phases of the N = n = 2 Renyi entropy 

In this appendix we argue that for N = n = 2, the only phase transition in the Renyi 
entroy is at x = |, using crossing symmetry as well as some properties of / which 
are unproven but checked numerically. (This was derived by a different method in 
[H].) Thus the Renyi entropy S A is always dominated by the vacuum block in one 
channel. For this calculation we will also assume that the OPE coefficients of light 
operators with A < A sap do not grow exponentially with c. 

The goal is to show that the contribution of heavy operators in the t-channel is 
exactly equal, to leading order in 1/c, to the vacuum contribution in the s-channel, 
and vice-versa. First we observe that the sum over light states is always dominated 
by the vacuum state at any value of x < |, not just near x = 0: 

Ope-W.A"*) ~ e-W°(/^) 0<x<±, (B.l) 

p|<5 p <l/24 

where 

Pn = -, 6 p =^. (B.2) 
nc nc 



26 



This follows from the fact that for n < i and any fixed x < |, f(5 p ,/3 n ,x) is an 
increasing function of <5 P in the range < 5 P < 1/24. This is straightforward to check 
numerically, and holds for any n. 

Now, crossing symmetry is the statement that the expansions in the two channels 
are equal: 

£ a p [ e -Wf(^M _ e -f = o . (B.3) 

p 

For each term in the sum, we want to keep only the dominant exponential. For this 
we need the behavior of the function 



D{6, f3, x) = f(S, 0, 1 - x) - f(5, p, x) . (B.4) 



When D is positive, the s-channel term in (B.3) dominates, and when it is negative, 
the t-channel term dominates. Numerically, we find the following: Depending on 
5, (3, there are three different behaviors for D as a function of x. D is monotonically 
increasing as a function of x for large 5, monotonically decreasing for small 5, and 
non-monotonic (with multiple roots) in some intermediate range. When n < 2, the 
intermediate range of 5 lies below 5 = 1/24, and when n > 2 it lies above 5 = 1/24. 
For n — 2, the intermediate range disappears; D is decreasing for 5 < 1/24 and 
increasing for 5 > 1/24. We caution that these statements rely on a numerical 
calculation of the conformal block near the region where the series diverges, x — > 1. 
This means that we have actually checked them only on the interval x G [e, 1 — e] 
with e ~ 0.01. We will assume that the statements hold exactly, but the numerical 
evidence is not sufficient to be entirely convincing. 



Returning to (B.3), if we take x < \ and n < 2, then these properties of D imply 



that the crossing equation becomes 

e -fM6 P ,M_ ^ a p e" W'/W-*) = , ( B .5) 

p|,5 p >l/24 

where we have dropped an irrelevant prefactor on the first term. This is exactly 
what we set out to show: the vacuum contribution in the s-channel is equal to the 
contribution of all heavy operators in the t-channel. Although this argument formally 
holds for 1 < n < 2, it tells us nothing about the entanglement entropy n — 1 because 
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we should not trust the analytic continuation from the single integer point n = 2. 
This argument fails for n > 2 because D is not monotonic in the range 5 > 1/24, 



so the dominant contributions to (B.3) do not separate cleanly into light operators in 
the s-channel and heavy operators in the t-channel. In fact, using the Cardy growth 
of states above 5 = 1/24, it can be shown that to avoid other phases for n > 2, a p 
must be exponentially suppressed in c around 5 ~ 1/24. It would be interesting to 
explore whether a p can be bounded some other way to demonstrate this suppression, 
perhaps along the lines of [50J. 
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